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In [Phys. Rev. A 58, 1833 (1998)] a family of polynomial invariants which separate the orbits 
of multi-qubit density operators p under the action of the local unitary group was presented. We 
consider this family of invariants for the class of those p which are the projection operators describing 
stabilizer codes and give a complete translation of these invariants into the binary framework in which 
stabilizer codes are usually described. Such an investigation of local invariants of quantum codes is 
of natural importance in quantum coding theory, since locally equivalent codes have the same error- 
correcting capabilities and local invariants are powerful tools to explore their structure. Moreover, 
the present result is relevant in the context of multipartite entanglement and the development of 
the measurement-based model of quantum computation known as the one-way quantum computer. 
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I. INTRODUCTION 

The theory of quantum error-correcting codes consti- 
tutes a vital ingredient in the realization of quantum com- 
puting, as these codes protect the vulnerable information 
stored in a quantum computer from the destructive ef- 
fects of decoherence. The most widely known class of 
error-correcting codes is that of the stabilizer codes, stud- 
ied extensively in e.g. 0,0,0 ■ An n-qubit stabilizer code 
is defined as a simultaneous eigenspace of a set of com- 
muting observables in the Pauli group, where the latter 
consists of all n-fold tensor products of the Pauli matrices 
and the identity. Equivalently, a code is described by the 
projector operator on this eigenspace. In the characteri- 
zation of the the error-correcting capabilities of quantum 
codes, two equivalence relations arise naturally on the set 
of corresponding projectors: two n-qubit quantum codes 
described by projectors p and p' are called globally equiv- 
alent, or just equivalent, if there exists a local unitary op- 
erator U E U(2)® n such that U pW is equal to p' modulo 
a permutation of the n qubits. If UpW — p', without 
any additional permutation, the codes are called locally 
equivalent. As globally equivalent codes have exactly the 
same error-correcting capabilities and vice versa, global 
equivalence is in fact the true equivalence of quantum 
codes. However, the structure of local equivalence is more 
transparent and insight in this matter already provides a 
lot of information about the structure of quantum codes. 
Therefore, much of the relevant literature tackles local 
equivalence and we will do the same in the following. 

This paper is concerned with the characterization of 
the local equivalence class of a stabilizer code p by means 
of local invariants. These are complex functions F(p) 
which remain invariant under the action of all local uni- 
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tary transformations, i.e., 

F{p) = F(UpU^) 

for every U £ U (2)®". In studying invariants, the general 
goal is to look for a minimal set of invariants which char- 
acterizes the local equivalence class of any given code. 
To obtain such a minimal complete set, it is well known 
4] that it is sufficient to consider functions F which are 
polynomials in the entries of p. These polynomial invari- 
ants form an algebra over C, as linear combinations and 
products of invariants remain invariants. Interestingly, 
the invariant algebra of {/(2)® n is finitely generated [5j 
and therefore the existence of a finite complete set of 
polynomial invariants is guaranteed. Although the prob- 
lem of pinpointing such a finite set is to date unanswered, 
progress has been made in the past years in construct- 
ing complete though infinite families of invariants. A 
natural approach is to consider homogeneous invariants, 
since any invariant can be written as a sum of its ho- 
mogeneous components, each of which needs to be an 
invariant as well. As the set of homogeneous invariants 
of fixed degree has the structure of a vector space, one 
wishes to construct a basis of this vector space degree 
per degree in order to obtain a generating (yet infinite) 
set of the invariant algebra. Although this goal has not 
yet been achieved, Grassl et al. |(| have constructed gen- 
erating (though nonminimal) sets of these spaces, using 
earlier work of Rains |7( • Their study of local invariants is 
general in the sense that is does not merely regard (pro- 
jectors associated with) quantum codes, but in fact ar- 
bitrary n-qubit density operators p. Closer inspection of 
their basic invariants when dealing with stabilizer codes 
is certainly appropriate, given the very specific structure 
of these codes and their associated projectors. Indeed, 
the stabilizer formalism has an equivalent formulation in 
terms of algebra over GF(2) and in this framework any 
stabilizer code of length n and dimension k is essentially 
described by an In x k binary matrix, called the gener- 
ator matrix of the code. It is therefore natural to ask 
what the structure of the basic invariants is in relation 
to this binary description. In the present paper we re- 



solve this issue: it is shown that the invariants of Grassl 
et al. are in a one-to-one correspondence with the di- 
mensions of certain linear subspaces over GF(2) which 
depend solely on the binary matrix description of a code. 
Hence, a complete translation of the invariants into the 
binary stabilizer framework is obtained. 

We wish to point out that the relevance of this inves- 
tigation stretches beyond the domain of quantum cod- 
ing theory: stabilizer codes with rank one projectors 
correspond to the class of pure states generally known 
as stabilizer states. The problem of recognizing local 
unitary equivalence of stabilizer states is of importance 
in the study of multipartite entanglement |a, li| and 
in the development of the one-way quantum computer, 
a measurement-based model of quantum computation 
which uses a stabilizer state as a universal resource 101. 



II. NOTATIONS 

First, we fix some basic notations which will be used 
throughout this paper. Seeing that stabilizer codes have 
descriptions both as projectors on a complex Hilbert 
space and as binary linear spaces, we will be dealing 
with algebra over the fields C and F2 := GF(2), where 
the latter is the finite field of two elements (0 and 1), 
where arithmetics are performed modulo 2. The set of 
p x q matrices over a field F G {C, F2} will be denoted by 
M pxq (¥), where p, q G Nq. To shorten notations, the set 
of square p x p matrices is denoted by M P (F). 

The group S r is the symmetric group of order r. For 
any n G No, 5" denotes the n-fold cartesian product 
of S r with itself, i.e., 5™ consists of all n-tuples II := 
(71-1, . . . , 7r„), where 7Ti G S r for every i = 1, . . . ,n. 



III. STABILIZER CODES AND LINEAR 
SPACES OVER GF(2) 

The Pauli group Q n on n qubits consists of all 4x4™ 
n-fold tensor products of the form a V\ <g> V2 <8> • • • <8> v n , 
where a G {±l,±i} is an overall phase factor and the 
2 x 2-matrices vi [i = 1, . . . , n) are either the identity <7o 
or one of the Pauli matrices 

a * = ( 1 ) ' °y = ( °i ) ' az = ( -1 ) • 

An n-qubit stabilizer S in the Pauli group is a subgroup 
of Q n which is generated by k < n commuting, indepen- 
dent and Hcrmitian observables Mj G Q n (i = 1, . . . , k). 
Here "independent" means that no product of the form 
Mf 1 . . .M£ k , where af» G {0, 1}, yields the identity ex- 
cept when all equal to zero. The stabilizer code 
associated with S is the joint eigenspace belonging to 



eigenvalue one of the k operators Mj. The numbers n 
and k are called the length and the dimension of the 
code, respectively. There is a one-to-one correspondence 
between the code associated with S and the matrix 

Mes 

as this operator is (up to a multiplicative constant) the 
projection operator which projects on the code space. 
The normalization is chosen such as to yield Trips) = 1. 

We now briefly discuss the binary representation of the 
stabilizer formalism (for literature on this subject, see e.g. 
PI 0| ) ■ Employing the mapping 

Co = 0-00 i-» (0, 0) 

cr x = 001 i-» (0, 1) 

v z = trio »-> (1,0) 

(T y = an ^ (1,1), (2) 

the elements of Q n can be represented as 2n-dimensional 
binary vectors as follows: 

<Ju lVl (8 • • • (8) cr«„„„ = (T(u,v) ^ (u, v) G Ff 1 , 

where (u, v) — (tti, . . . , u„, V\, . . . , v„). This parameteri- 
zation establishes a group homomorphism between Q n , ■ 
and F|™, + (which disregards the overall phases of Pauli 
operators). In this binary representation, two Pauli op- 
erators a a and Cb, where a, b G Ff", commute if and only 
if a T Pb = 0, where the 2n x 2n matrix 



defines a symplectic inner product on F|". Therefore, 
a code of length n and dimension k corresponds to 
a ^-dimensional linear subspace of F| n which is self- 
orthogonal with respect to this symplectic inner product, 
i.e., a T Pb = for every a, b in this subspace. Given a 
set of generators of the stabilizer, we assemble their bi- 
nary representations as the columns of a full rank 2n x k 
matrix 5, which is referred to as a generator matrix of 
the stabilizer subspace. This generator matrix satisfies 
S T PS = from the symplectic self-orthogonality prop- 
erty. The entire binary stabilizer subspace (or code space) 
Cs consists of all linear combinations of the columns of 
S, i.e., it is equal to 

C s := {Sx I x G F2} . (3) 

It is important to notice what happens when some of 
the qubits in the system p$ are traced out. For every 
u> C {1, . . . , n}, the partial trace Tr w c p$ =: p% yields a 
stabilizer code on \u>\ qubits, where uj c is the complement 
of u> in {1, ... ,n}. Using the definition it follows that 



As all three Pauli matrices have zero trace, the sum can 
be taken over all M G S which are equal to the identity 
(To on the ith tensor factor for every i G uj c . Defining 
the support supp(M) of any M = a V\ ® • • • <g> v n 6 Q n 
by the subset of those i € {l,...,n} such that V{ ^ 
(To, the sum in (0J runs over the subgroup of all M G 
S such that supp(M) C u>. For such M, the partial 
trace Tr w c M removes the tensor factor ao on positions 
i G lo c . Transferring the definition of supp to the binary 
representation of Q n , the binary code space Cg of pg is 
obtained by considering the subspace of those y G C$ 
such that supp(y) C cj and removing form these y's the 
components (yi,y n +i) — (0,0) for every i G oj c . 

Stabilizer states and graph states. If the dimension of 
a stabilizer code is equal to its length, i.e., if k = n, 
then the code is called self-dual. It is easy to see that 
the code space of a self-dual code is one-dimensional or, 
equivalently, p$ — for some pure state The 

states \tjj) obtained in this way are known in the literature 
as stabilizer states. By definition, a stabilizer state is the 
unique simultaneous eigenvector with eigenvalue 1 of a 
set of n commuting and independent Pauli operators. A 
subset of the class of stabilizer states which will be of 
particular interest in our investigation is constituted by 
the so-called graph states pToL fll| . For these states, the 
defining eigenvalue equations can be constructed on the 
basis of a graph: when G is a simple graph on n vertices 
with adjacency matrix 9 one defines n (commuting) 
Pauli operators 



Ki 



71 



fc=i 



where o4 , o-\ l) > o^ 1 are the Pauli operators which have 
resp. o~ x ,o~y, o~ z on the ith position in the tensor product 
and the identity elsewhere. The graph state \G) is the 
stabilizer state associated with the operators Kj, j = 
1, . . . ,n. The rank one projector \G)(G\ is denoted by 
pa- Note that the binary code space of a graph state 
\G), for a graph G with adjacency matrix 6, is generated 
by 



Local Clifford operations. The Clifford group C\ on one 
qubit is the normalizer of Q\ in U(2), i.e. it is the sub- 
group of 2 x 2 unitary operators which map Q\ to itself 
under conjugation. The local Clifford group C l n := Cf n 
on n qubits is the n-fold tensor product of C\ with itself. 
When disregarding the overall phases of the elements in 
Qi, it is easy to see there exists a one-to-one correspon- 
dence between the one-qubit Clifford operations and the 
6 possible invertiblc linear transformations of F|, since 
each one-qubit Clifford operator performs one of the 6 
possible permutations of the Pauli matrices and leaves 
the identity fixed. Generalizing to n-qubit local Clifford 



operations, it follows that each U G C l n corresponds to a 
matrix Q G M2 n {^2) of the block form 



Q 



A B 
C D 



where the n x n matrices A,B,C,D are diagonal. We 
denote the diagonal entries of A, B, C, D, respectively, 
by <2j, bi, Ci, di, respectively. The n submatrices 



Q 



M — 



bi 
di 



correspond to the tensor factors of U. It follows from the 
above discussion that each of the matrices is invert- 
ible. We denote the group of all such Q by C n . It follows 
that two stabilizer codes ps, ps> with generator matrices 
S, S', respectively, are equivalent under the local Clifford 
group if and only if there exists Q G C l n such that 



Cqs = Cs> 



(5) 



To see this, simply note that UpslJ 1 = ps> for some 
U G C' n if and only USW =5'. 

Finally, for our investigation it is important to note 
that any stabilizer state is locally equivalent to a (gen- 
erally nonunique) graph state under the local Clifford 
group 0. 



IV. INVARIANTS, PERMUTATIONS AND 
BINARY TREES 



In this section, we recall the constructions of basic 
polynomial invariants reported in refs. 0, 0j- 

Let p G (C) be an n-qubit density operator. Any 
homogeneous polynomial F(p) of degree r in the entries 
of p can be written as a trace 

F(p) = Tr (A F ■ p^) 

for some matrix Ap G M2»«-(C). To see this, simply 
note that the tensor product p® r contains all monomials 
of degree r in the entries pij . The coefficients of these 
monomials in the polynomial F are encoded in the entries 
of Af. Consequently, F(p) is an invariant of U(2)® n if 
and only if there exists an Ap such that 



[Af,^ 







(6) 



for every U G J7(2)® n . Therefore, the study of invari- 
ant homogeneous polynomials of fixed degree r is trans- 
formed to the study of the algebra of matrices A p which 
satisfy (©■ It was shown by Rains that a set of matri- 
ces which linearly generate this algebra can be obtained 
in a one-to-one correspondence with the group 5™ as fol- 
lows: let II = (//, v, £, . . . ) G <S™ be an n-tuple of per- 
mutations. The matrix Tu G M2»r(C) is defined as the 
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FIG. 1: Binary tree on 10 nodes with maximal right paths 
(1, 3, 9, 10), (2), (4, 7, 8) and (5, 6). Note the canonical way in 
which the nodes are labelled. 



permutation matrix which acts on (C 2 )® r by permuting 
the r copies of the ith qubit according to the zth permu- 
tation of II, i.e., Xn maps a tensor 



tree B is defined by the product of cycles 



"0H ,jl,kl...\ l2,]2,k 2 



i r ,j r ,K... G (C 2 ")« 



to 



V"V(1) >Ml) —i V(2) >J»(2) < k Z 



If II ranges over all elements in 5™ , the matrices Tn lin- 
early generate the algebra defined by ©• Therefore, one 
obtains a generating set of basic invariants I r ji of degree 
r, where 



IrM :=Tr (T u -p® r ). 



(7) 



However, linear dependencies within the set of matrices 
Tn do exist. In ref. |g Grassl et al. improved the above 
result, as the authors presented a method which is able 
to pinpoint within the set {T,n}n a linearly independent 
subset for every r. Their approach was to consider binary 
trees and to associate with every binary tree B on r nodes 
a permutation tt(B) € S r . Enumeration of all possible 
n-tuples of permutations obtained in this way then yields 
a basis of the space of matrices Ap which satisfy ©. We 
now repeat the details of this construction. 

A (labelled, ordered and connected) binary tree B on 
r vertices is a special instance of a simple, oriented and 
connected graph, i.e. it consists of a set of vertices or 
nodes V = {1, . . . , r} which can be connected by arrows 
according to a number of prescriptions. If there is an 
arrow from a node / G V to a node s G V then / is 
called the father of s and, conversely, s is a son of /. In 
a binary tree, all nodes but one have exactly one father. 
The one node without father is called the root of the tree. 
Furthermore, every node has at most two sons (called 
left and right son, respectively). The labelling of the 
r nodes is obtained by traversing the tree in the order 
root - left subtree - right subtree. A maximal right path 
p in a binary tree B is an ordered tuple of nodes p = 
(uq, vi, . . . , v s ) such that vq is not the right son of any 
node of B, Vi is the right son of for i = 1, . . . , s and 
v s has no right son. An example of a labelled binary tree 
is given in Fig. 1. 

Denoting by 71(B) the set of all maximal right paths 
of B, the permutation ir(B) associated with the binary 



n(B) 



n 



(v Vi . ..v s ). 



e K(B) 



Note that n(B) G S r whenever B has exactly r nodes 
and that there is a one-to-one correspondence between 
B and tt(B). The set of all permutations obtained in 
this way is denoted by V r . According to the result in 
Kg, the invariants {I r ,n}, where IT = (m,..., n r ) G V" 
ranges over all n-tuples of permutations in V r , forms a 
basic set of homogeneous invariants of degree r. 

To conclude this section, we state some definitions re- 
garding binary trees, which will be used below. Let B 
be a binary tree on r nodes. The start st(p) of a path 
p = (vq, vi, . . . , v s ) € Tt{B) is the element vq and the fin- 
ish fin(p) is the element v s . The length of p is the num- 
ber s + 1 . By an expression of the form "i 6 p" is meant 
that the node i belongs to the set {vq,vi, ... ,v s } (note 
that, due to the canonical labelling of the nodes, there 
is a one-to-one correspondence between p and the set 
{vo, vi, . . . , v s }). For every node i, the path p(i) e 71(B) 
is the unique maximal right path such that i € p(i). 

Let B have t := \1Z(B)\ maximal right paths p±, . . . ,pt, 
which we suppose to be ordered in such a way that st(pi) 
< st(j>2) < ■ ■ ■ < st(pt). The columns {Rb)j of the matrix 
Rb G Af rxt (F 2 ) are defined by: 



(Rb)j = ^ ej, 



(8) 



for every j G {l,...,i}, where e; is the ith canonical 
base vector in F£. The columns (Db)j of the matrix 
Db G M r (F 2 ) are defined by: 



{Db)j 



E « 



(9) 



for every j G {1, ...,r}. Finally, the linear space Vb 
consists of all x G F£ such that J2ie P Xi = for every 
p G 71(B) (i.e., Vb is the null space of the matrix -R^). 



MAIN RESULT AND DISCUSSION 



We are now in a position to state the central result of 
this paper. 

Theorem 1 Let p$ be a stabilizer code of length 
n and dimension k with generator matrix S. Let Sf 
(i = 1, . . . ,n) be the 2 x k submatrix of S obtained by 
selecting the ith and the (n + i)th row of S. Fix r G No, 
let Bi, B2, ... ,B n be n binary trees on r nodes and let 
IT G V™ be the associated n-tuple of permutations. Then 



logi I r ,n(ps) ~ dimw 2 ker 



R B 2 



R B„ 



'ST 



(10) 



where ~ denotes equality up to an additive constant in- 
dependent of ps ■ 

Theorem 1 shows how the information contained in 
the invariants I r ji can be recuperated within the binary 
representation of the stabilizer formalism. The fact that 
a translation into the binary framework is possible is of 
course not unexpected, as a stabilizer code is, up to in- 
formation about the overall phases of its stabilizer ele- 
ments, defined by its generator matrix. Moreover, these 
phases do not play a role in determining the (local) equiv- 
alence class of a code 0. However, the simple form of 
the result (|10|l is remarkable: an invariant I r: u is in a 
one-to-one correspondence with the dimension of a bi- 
nary linear space which depends only on the generator 
matrix S - and this in a very transparant way. Addition- 
ally, it is interesting to notice the explicit way in which 
the n-tuple of binary trees appear in the result: every 
matrix corresponding to the ith binary tree, is cou- 
pled via a tensor product to the matrix Si, which is the 
subblock of S containing the information about the ith 
qubit. Finally, we note that the r.h.s. of (|10|l can be 
computed efficiently via a calculation of the rank over F2 
of the matrix 

' RT B,® S l' 

Rg 2 (g> #2 ^ 

Before proving theorem 1 in section VI, we investigate 
the invariants <|1(J|I in more detail. We start with the in- 
variants of smallest nontrivial degree, i.e. r = 2. There 
are exactly two binary trees on 2 vertices, as node 2 can 
either be the right or the left son of node 1. Equiva- 
lently, there are two possible matrices Rb according to 
definition JSJ, namely 

I ?] and [11] T , (12) 

where the identity matrix corresponds to the tree where 2 
is the left son of 1. Now, consider an n-tuple (Bi, . . . , B n ) 
of binary trees and the corresponding n-tuple of permu- 
tations II G V2 ■ Let uj C {1, . . . , n} denote the set of 
all i such that = [1 1] T - note that every n-tuple of 
binary trees on 2 nodes corresponds uniquely to such a 
set uj. Using the notation Si as in theorem 1, this implies 
that 

R T Bi ® Sj = [ST Sf] 
whenever i£w and 

otherwise. Therefore, the null space of the matrix 
consists of all vectors (x,x') G F| fe such that 

Si(x + x') = for every i G uj 
SjX = Sjx' = for every j G uj c , (13) 



where uj c is the complement of uj in {l,...,n}. Note 
that l|13|) implies that Sx = Sx' and therefore x = x' , 
since S has full rank. Thus, the solutions of are in 
a one-to-one correspondence with the linear subspace of 
F§ of those x satisfying SjX = for every j G ui c . The 
linear mapping 0s : Fij — > F|™ defined by the matrix S 
maps the space of such x's to the space of vectors y = Sx 
which satisfy yj — y n +j — for every j G uj c . As S has 
full rank, the mapping <fis is injective and the spaces of 
the x's and the y's have equal dimension. Recalling that 
the support supp(w) of any v G F2™ is the subset of those 
i G {1, . . . ,n} such that {vi,v n+ i) ^ (0,0), we can state 
that the supports of the y's lie within the set uj. Thus, 
we have shown that 

log 2 h,n ~ dim {y G C s \ supp(y) C uj}. (14) 

This is clearly a more insightful presentation of the invari- 
ants 22, n than IjlUII . as (|14f) relates invariants to the di- 
mensions of the subspaces Cg of the code space Cs- More- 
over, a similar argument as above can be made to obtain 
an analogous presentation of the invariants of higher de- 
gree: 

Theorem 2 Let S be a 2n x k generator matrix of a 
stabilizer code and let the subblocks Sf be defined as in 
theorem 1. Fix r G No and let B\, . . . , B n be n binary 
trees on r vertices. Let 1Z = TZ(Bi) U • • • U lZ(B n ) denote 
the set of all maximal right paths of these trees. For every 
p G 1Z, let uj p denote the subset of all i G {1, . . . , n} such 
that p ^ TZ(Bi). Then the dimension of the kernel of the 
matrix 1) 11)) is equal to the dimension of the space 

{(y« G C s x ... x C s \ 

supp C^^y^) Q <-J p , for every p G 72.}. 

(15) 

The proof of theorem 2 is omitted, as it is a straightfor- 
ward generalization of the considerations made above for 
the invariants of degree 2. A number of properties of the 
invariants immediately follow from theorem 2: e.g., con- 
sider an invariant / r +i,n such that its n binary trees Bi 
have got a common maximal right path po of exactly one 
element i , i.e. po = (io) £ T^(Bi) fl ■ • • D 7Z(B n ). Con- 
sidering i|15|) for this invariant, it follows that uj Po = 0. 
Consequently, supp(?/ 10 )) C and therefore y^ = for 
every y^ in (|15fl . Thus, the invariant / r +i,n is equal to 
an invariant of degree r corresponding to binary trees Bi 
which are obtained by deleting the node io from the orig- 
inal trees Bi. Reversing the above argument shows that 
any invariant of degree r can be written as an invariant 
of degree r + 1, thereby showing that 

{2V,n}n C {ir+i,ri'}n' 

for every r, where II (II') ranges over all elements in T 5 " 
(Pr+i)- Moreover, this argument can be generalized to 
show that an invariant corresponding to a tuple of binary 
trees which have some subtree in common, can be written 
as the product of invariants of lower degree. 



VI. PROOF OF THEOREM 1 



The proof of theorem 1 is given in two main parts. 
We show in subsection A that it is sufficient to prove 
the theorem for graph states. The proof of theorem 1 
for graph states is subsequently given in subsection B. 
For the remainder of this section, we fix r and consider 
an n-tuplc II = (jri, . . . , ir n ) of permutations 7Tj € V r 
and n binary trees Bi on r vertices corresponding to the 
permutations 7^. We denote B := {B\, . . . , B n ). 



A. Reduction to graph states 



Suppose that theorem 1 holds for all graph states. Let 
ps be a stabilizer code of length n and dimension k with 
generator matrix S. For large enough m > n, there exists 
a stabilizer state \tf>) on m qubits such that ps can be 
obtained from \tp) by tracing out the qubits n + 1, . . . , m, 
i.e. 

ps = k/>>(V| {1 '"- n} 

using the notation of section 3. Furthermore, |^) is equiv- 
alent to some graph state \G) under the local Clifford 
group. Denoting co = it follows that p$ is 

locally equivalent to Pq. Letting S' be the generator ma- 
trix of |G), this last fact translates into the binary picture 
as 



C QS = Cg, 



(16) 



for some Q e C l n . Now, let IT G V™ be the m-tuple of 
permutations which is obtained by appending to II m — n 
times the identity permutation (which belongs to V r ) and 
let B' = (Bx, . . . , B n , B , . . . , B ) be the associated m- 
tuple of binary trees; here Bg is the binary tree with r 
maximal right paths (i), corresponding to the identity 
permutation. The crucial observation is now 

Ir.n(Ps) = Ir,w{4>)- 

This identity can easily be verified by using the definition 
of I r ji>. Moreover, I r jii(ip) = Ir.n'ipc) since the states 
\ip) and \G) are locally equivalent. We can now apply 
theorem 1 and find 



log 2 -fr,n(ps) = log 2 / r! n'(pG) 





" R T Bi i 




ker 


R T Ba % 






- Rl f < 


5 ^rJ - 



(17) 



Applying theorem 2 to the generator matrix S' and the 
binary trees £>', H17|) is equal to the dimension of 



{(y 



1| ) ...,y M ) G C s < x ••• x C s ,\ 

supp \ / j y^') C Lu p , for every p G 7Z'} 



(18) 

where W = K(B X ) U ■ ■ • U K(B n ) U K(B ). As the last 
m — n trees in the m-tuple B' are equal to -Bo, every yV> 
in l|18|) has supp(y(^) C uj. Therefore, the dimension of 
(fTHfl is equal to the dimension of 



{(■ 



supp (y j x^') C cjp, for every p G TZ} 



(19) 



where now 1Z = TZ{B{) U • ■ • U lZ(B n ). Finally, we recall 
the identity (|16|l and note that l|19l) remains invariant if 
Cqs is replaced by Cs- A last application of theorem 2 
yields 



log 2 Ir,n(As) ~ dim F2 ker 



which is the desired result. 



R B 2 



si 
si 



Rl ®SI 



(20) 



B. Proof of theorem 1 for graph states 

Fix a graph G on n vertices with adjacency matrix 9 
and the generator matrix 

3 = 

which has 2x n subblocks Sf defined as in theorem 1 by 

sT = 



The proof of theorem 1 for the case where ps = pa is 
structured as follows: in lemma 3 we show that the in- 
variant I r ,Yi{pG) is equal to a sum of the form 



(21) 



xev 



where N is a normalization factor independent of G, V is 
a linear subspace of M„ xr (F2) and Q is a quadratic form 
on M„ xr (F2). Preliminary material used to prove this 
result will be gathered in lemmas 1 and 2. In lemma 4, 
we subsequently show that the form Q is in fact identical 
zero on the space V, which implies that the sum in l|21|) 



is (up to the normalization) equal to the cardinality of 
V. Finally, this cardinality is related to the r.h.s of 110(1 
and the proof of theorem 1 is completed. 

It will be convenient to work with a real variant of the 
set of Pauli matrices (as in ref. 13]), defined by 



Too 
roi 
no 

m 



ia n = 



1 

-1 



(22) 



which we will call the tau matrices. Analogous to the 
notation introduced in section 3, n-fold products of tau 
matrices are represented as 

where (u,v) = (ui, ■ ■ ■ ,u n , vt, . . . , v n ) £ F|™. We now 
prove a useful parameterization of the projector pa- 

Lemma 1 The projector pa can be parameterized as 
follows: 



PG 



on ^ ; ( ^ 



kg(x) 



(6x,x) : 



(23) 



xeF" 



where kg is the quadratic form over F 2 associated with 9, 
i.e., k s (x) = J2i< 3 S lj x l x J . 

Proof: the state \G) is defined by the n relations 
T(g j e ,.)|G) = |G), where 9j is the jth column of 9 and 
ej is the jth canonical basis vector of F 2 . The stabilizer 
of \G) consists of all products 

n 

i=i 

where x — (xi, . . . , x n ) £ FJj. After a repeated applica- 
tion of the multiplication rule |l3j 

/ -1 \u T u' 

'(u,v)'(u',v') I ^} ' (u+u' ,v-\-v') ) 

where u, u', v, v' £ FJ, we arrive at 

Mx = (_l)M*) T(ex :c) , 

Since pa — ^ Sksf™ -^as) we obtain the result. □ 

Lemma 1 will be used below to compute the invariant 
Ir,n(pG)- After plugging l(2*3l in J7J), we will be dealing 
with expressions of the form 



Tr (T n n ® r 2 



Tr), 



(24) 



where the Tj's are themselves n-fold tensor products of 
the tau matrices, i.e. £ C/„ for every i = 1, ...,r. 



A closer look at these expressions beforehand is appro- 
priate. To this end, let X,Y, . . . be any r operators in 
M 2 (C)® n , i.e., 

X = Xi ® • ■ • ®X n , 
Y = Yi <8>---®y n , 

where Xi,Yi, ■ ■ ■ £ Af 2 (C). Furthermore, for any tt £ V r 
we denote 

A-,tt(^j V, ■ ■ ■ ) := ^2 ^»2V(2) • • • > 



where [/, V, . . . are r arbitrary 2x2 matrices. Using the 
definition of Tn, it is then easy to check that 

n 

Tr {T R X®Y®...) = l[A rtWi (X i ,Y i ,...). (25) 

i=l 

It follows that (|24|l is a product of n factors of the form 

At^-k {t~u\ V21 ' ' ' 1 7~u r v r ) — ■ A rlx (li, u) , 

where u = (ui, . . . ,u r ) and i> = (v\, . . . ,v r ) € F 2 . Ex- 
pressions of this type are calculated in lemma 2: 

Lemma 2 Let tt eV r be a permutation corresponding 
to a binary tree B and let (u,v) £ F 2 r . Let the matrix 
Db and the space Vb be defined as in section 4- Then 



'tr—dim Vf 



(_i)« u B v i f u ^ ve Vb 

otherwise 



Proof: First, note that the entries of the 1-qubit oper- 
ators Tab can be parameterized as 

(r ab ) x . v = (-l) a ^S x+Vlb , 

where a, 6, x, y £ F 2 . Using this formula in the definition 
of A Ti „, we obtain 

Ar^iu, fc>) — ^ y ( 1) ^ ^xi+Kjrfij ,vi ■ ■ ■ 0~x r -\-x^( r ) ,v r • 

xefi; 



Equivalcntly, the sum runs over all x £ F 2 which 
lie in the affine subspace determined by the equations 
Xi + x^u\ = Vi for all i — 1, . . . , r. However, this system 
of equations may not be consistent: indeed, one can eas- 
ily show that a solution exists iff v £ Vb. Whenever this 
is the case, the solutions are given by x = xq + x' , where 
xq = D^v + v and x' satisfies x\ -I- x' n u\ = for every 
i = 1, . . . , r. Moreover, the space of all such x' is the or- 
thogonal complement of Vb with respect to the standard 
inner product in F 2 (or equivalently, the column space of 
Rb), as one can verify. Therefore, we have 



A r ^{u,v) 



Y.x>ev^- l ) uT{DTBV+x ' ] if«eU B 



otherwise 



Furthermore, the sum X^x'eV^-( — -*-)" x 1S ec l ua l to 
2r-dim v B jf u g y B anc [ zer0 otherwise. This proves 
the result. □ 

We now proceed in calculating the invariant I r j\{po). 
Using lemma 1 in (J7J), we find that I r ,n(pG) is equal to 
the sum 

J_ y ((-ijEUM^ix 

Tr (T n T (ex (l) ia .(l)) O ■ ■ ■ <g> T (fe (r) i;c (r)))} 

(26) 

First, denoting by 'Plows (#) the strictly lower triangular 
part of and writing 

X := [arW|...|x«] G M nxr (F 2 ), 

we obtain the shorthand notation 

r 

^Mz w ) = TrX T P lows (#)X 



Secondly, the trace in (|26|) splits into a product of n fac- 
tors as in (|25ll . each of which can be calculated by employ- 
ing lemma 2. The calculation is straightforward. Defin- 
ing for every X G M nxr (F2) the matrix Xg G M„ xr (F2) 
by 



fe=i 



one finds that 1)26(1 is equal (up to a normalization inde- 
pendent of G) to 

^(_!)TrX r P lows (e)X+TrXjeX ; ^ 

where the sum runs over all X such that 



S : 







(28) 



for every i G {1, . . . ,n} and p G IZ(Bi). We will denote 
the space of all such X by Vg(G). We have proven: 

Lemma 3 The invariant Irjiipc) can be written as 

(29) 

where N is a normalization factor independent of G and 
the definitions o/Vg(G) and Xg are as above. 

The last part of our argument consists of showing 
that the quadratic form Q(X) := Tr X T V\ ovls {9)X + 
Tr Xg OX is zero on the space Vg(G). Once this 
result is shown, the proof of theorem 1 is immedi- 
ate: indeed, if Q(X) = for every X G Vb{G) then 



log 2 ^r,n(pG) 
X = [x<U\...\x<r 
tors X = (x^, . . . , x^ ) S FY which are exactly the 
elements in the null space of the matrix 



dim Vg(G). Moreover, the matrices 
'] G Vg(G) can be reshaped as vec- 



R T B2 ®S? 



•si 



R B„ 



(30) 



Clearly, the spaces of the X's and the X's have the same 
dimension and the proof of theorem 1 is thus completed. 
We now show that Q = on the space Vg(G): 

Lemma 4 Q(X) = for every X G V B (G). 

Proof: Let X = [xW \... \x^} be an element of V B (G). 
Recall that by definition (|28|l this entails that 



■ U) 



(31) 



for every i G {1, . . . , n} and p G IZ(Bi). In particular, 

for every i G {1, . . . , n} and for every p G IZ(Bi), where 
a;0) = (xi \ . . . ,Xn)- Consequently, 



(33) 



5>«> = o. 

Now, consider the first term of Q(X): 
Qi := Tr X T ^ lows (0)X 

r 

= J2x^ T V lows (6) x®. 

3 = 1 

Substituting — Y^=i x ^ (from (|33|l 'l. an easy cal- 
culation shows that 



(34) 



r-1 



.(AO 



,0') 



Let C { 1 , 



i} consist of all k G {1, . . . ,j — 1} which 



belong to a maximal right path p of such that fin(p) > 
j. Then, denoting := St=i l|32() implies that 



(.;) 



„(fc) 



The second term of Q{X) is 

Q 2 



Tr OX 



(35) 



where is the jth column of Let rjij consist of 
all k G {1, . . . , j — 1} which belong to the unique maximal 
right path of Bi which contains j. It then follows from 
the definition of that z^ = J2kevij^{j} x i ■ Note 
that rj ir U {r} G TZ(Bi) and therefore z{ r ' = for every 
i = l,...,n from (32). Thus, z (r) 
above results, we obtain 



0. Combining the 



r-1 



Qi + Q 2 = E^ (i) + zO) ) T0 ( 36 ) 



where 



# + * W = 




xf. (37) 



In the sum (|3T|> - every A: G Wjj n gives rise to a double 

appearance of the term a;^ and consequently all such 
terms vanish. As rjij C wy, we obtain 



,00 



z9 ] = 



E 



,,0) 



(38) 



When using (|!|SJ) to calculate the terms x\ J ' in (jSSJ) 
do not contribute to the sum, as they give rise to terms 

in H36|) . which are equal to zero since 9 is 
symmetric. Thus, defining the vectors «w by 



E 



„(fc) 



,k£uiij\rn 



(|3l))) becomes 



.(i) J 



rO') 



(39) 



Note that the set \77.y consists of all k G {1, . . . ,j — 1} 
which belong to some path p ^ p(j) in IZ(Bi) such that 
fin(p) > j. We now show that whenever j and i belong 
to the same maximal right path of one has u)ij \ r]ij = 

^ii\ilu and consequently iif = uf . To see this, fix i and 
consider arbitrary nodes j and I which lie on the same 
maximal right path of Bi. Without loss of generality 
we can assume that j < I. Denoting by j' the right 
son of j, we prove that u)\j \ rjij — LOiy \ rjij''. indeed, 
if j' = j + 1 then j does not have a left son (due to 
the canonical labelling of the nodes) and the assertion 
follows trivially; if on the other hand j' > j + 1 then j 
has a left subtree. However, the maximal right paths p 
in this subtree do not contribute to u>if \ rjij/ , as they all 
satisfy fin(p) < j' (which is again due to the canonical 
labelling of the nodes). Therefore ajy\rte = w,,-/ \rjiji and 
iteration of this argument shows that \ rjij = oju \ rju . 



We will now use the above property of the u^'s to 
show that Q(X) = 0. Let us consider the first binary 
tree B\ and suppose that (1, 2, 3) is a maximal right path 
of this tree. This example is chosen for notational conve- 
nience, but the argument will work for any maximal right 
path in any tree. Thus, we have — = = u. 



Denoting = (u, v^ ) for j 
terms in (|39|l are 



1,2,3, the relevant 



= (u,v^) T 8 + (u,vW) T 6 + ( u ,v^) T e x® 
= (u, Of^+^+^l 
+ (0,v^) T 8 + (0,v^) T 9 x<® + (0,v^) T 8 x (3) 

In the r.h.s. of the last equality, the first term is equal to 

u0[(xM +x& +x ^)) i 

which is equal to zero from l|31|l . since (1,2,3) is a maxi- 
mal right path of B\. Applying this argument to all the 
maximal right paths of the trees in IB shows that indeed 
Q = on the space Vb(G). This ends the proof. □ 

VII. CONCLUSION 



In this paper, we have considered a complete family of 
local invariants of stabilizer codes and we have given a 
translation of these invariants into the binary representa- 
tion of the stabilizer formalism. In particular, we have re- 
lated invariants to dimensions of binary subspaces which 
depend only on the generator matrix of a code. The aim 
of this investigation is mainly to provide a tool to study 
the structure of equivalence classes of codes. We note 
that some important issues in the present matter remain 
to be settled: firstly, it is to date not clear how a finite 
complete set of invariants can be constructed, i.e., what 
the minimal degree r is such that the values of the invari- 
ants of degree smaller than r determine the local equiv- 
alence class of any stabilizer code. Secondly, there is the 
question whether it is sufficient to consider local Clifford 
operations in order to recognize local equivalence of sta- 
bilizer codes. In other words, are two stabilizer codes 
locally equivalent if and only if they are equivalent under 
the local Clifford group? We believe that the results in 
this paper are a significant step towards answering these 
questions. 
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